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Abstract 


The production of forward jets separated by a large rapidity gap at LHC, the so-called 
Mueller-Navelet jets, is a fundamental testfield for perturbative QCD in the high-energy 
limit. Several analyses have already provided with evidence about the compatibility of 
theoretical predictions, based on collinear factorization and BFKL resummation of en¬ 
ergy logarithms in the next-to-leading approximation, with the CMS experimental data 
at 7 TeV of center-of-mass energy. However, the question if the same data can be described 
also by fixed-order perturbative approaches has not been yet fully answered. In this paper 
we provide numerical evidence that the mere use of partially asymmetric cuts in the trans¬ 
verse momenta of the detected jets allows a clear separation between BFKL-resummed 
and fixed-order predictions in some observables related with the Mueller-Navelet jet pro¬ 
duction process. 
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Figure 1: Mueller-Navelet jet production process. 

1 Introduction 

It is widely believed now that the inclusive hadroproduction of two jets featuring transverse 
momenta of the same order and much larger than the typical hadronic masses and being 
separated by a large rapidity gap Y, the so-called Mueller-Navelet jets [T], is a fundamental 
testheld for perturbative QCD in the high-energy limit, the jet transverse momenta providing 
with the hard scales of the process. 

At the LHC energies, the theoretical description of this process lies at the crossing point of 
two distinct approaches: collinear factorization and BFKL [2] resummation. On one side, at 
leading twist the process can be seen as the hard scattering of two partons, each emitted by 
one of the colliding hadrons according to the appropriate parton distribution function (PDF), 
see Fig. Collinear factorization takes care to systematically resum the logarithms of the hard 
scale, through the standard DGLAP evolution |3] of the PDFs and the hxed-order radiative 
corrections to the parton scattering cross section. 

The other resummation mechanism at work, justihed by the large center-of-mass energy 
a/s available at LHC, is the BFKL resummation of energy logarithms, which are so large 
to compensate the small QCD coupling and must therefore be accounted for to all orders of 
perturbation. These energy logarithms are related with the emission of undetected partons 
between the two jets (the larges s, the larger the number of partons), which lead to a reduced 
azimuthal correlation between the two detected forward jets, in comparison to the hxed-order 
DGLAP calculation, where jets are emitted almost back-to-back. 

In the BFKL approach energy logarithms are systematically resummed in the leading log¬ 
arithmic approximation (LLA), which means all terms (Q;sln(s))”, and in the next-to-leading 
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logarithmic approximation (NLA), which means resummation of all terms ln(s))”. The 
process-independent part of such resummation is encoded in the BFKL Green’s function, obey¬ 
ing an iterative integral equation, whose kernel is known at the next-to-leading order (NLO) 
both for forward scattering (he. for t = 0 and color singlet in the t-channel) [U [S] and for any 
hxed (not growing with energy) momentum transfer t and any possible two-gluon color state 
in the t-channel mm 

To get the cross section for Mueller-Navelet jet production and other related observables, 
the BFKL Green’s function must be convoluted with two impact factors for the transition from 
the colliding parton to the forward jet (the so-called “jet vertices”). They were hrst calculated 
with NLO accuracy in [S] and the result was later conhrmed in uni. A simpler expression, 
more practical for numerical purposes, was obtained in m adopting the so-called “small-cone” 
approximation (SGA) [T^ IT3] . he. for small jet cone aperture in the rapidity-azimuthal angle 
plane. A critical comparison between the latter result and the exact jet vertex in the cases of 
kf and cone algorithms and their “small-cone” versions has been recently carried out in [H]. 
We stress that, within the NLO accuracy, the jet can be formed by either one or two particles 
and no more, so that the argument given in na about the non-infrared-safety of the all-order 
extension of the jet algorithm used to obtain the SGA jet vertex in [TT] does not apply here. 

The BFKL approach brings along some extra-sources of systematic uncertainties with re¬ 
spect to the hxed-order, DGLAP calculation. First of all, in addition to the renormalization 
and factorization scales, fipt and fip, which appear also in DGLAP, there is a third, artihcial 
normalization scale, usually called Sq; which must be suitably hxed. Moreover, there is com¬ 
pelling evidence that choosing for these scales the values dictated by the kinematics of the 
process is not necessarily the best choice when the BFKL resummation is at work. It is well 
known, indeed, that the NLO BFKL corrections for the n = 0 conformal spin are with opposite 
sign with respect to the leading order (LO) result and large in absolute value. This happens 
both to the NLO BFKL kernel and to the process-dependent NLO impact factors (see, e.g. 
Ref. [I5], for the case of the vector meson photoproduction). This calls for some optimization 
procedure, which can consist in (i) including some pieces of the (unknown) next-to-NLO correc¬ 
tions, such as those dictated by renormalization group, as in collinear improvement [I6], or by 
energy-momentum conservation HU, and/or (ii) suitably choosing the values of the energy and 
renormalization scales, which, though arbitrary within the NLO, can have a sizeable numerical 
impact through subleading terms. Gommon optimization methods are those inspired by the 
principle of minimum sensitivity (PMS) [18], the fast apparent convergence (FAG) [T9| and the 
Brodsky-LePage-Mackenzie method (BLM) |2nj . 

This variety of options reflects in the large number of numerical studies of the Mueller- 
Navelet jet production process at LHG, both at a center-of-mass energy of 14 TeV [211122| |23| 
and 7 TeV [241 124| [2^ 127] . All these studies were concerned with the behavior on Y of azimuthal 
angle correlations between the two measured jets, i.e. average values of cos(n0), where n is 
an integer and 0 is the angle in the azimuthal plane between the direction of one jet and the 
opposite direction of the other jet, and ratios of two such cosines [281 . particular, one of 
these analyses [25], based on the use of a collinearly-improved and energy scales optimized a la 
BLM, found a nice agreement with GMS data |2U] . 
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In a recent paper 123, some of us stressed that another important source of systematics 
on should be aware of is the “representation uncertainty”, deriving from the freedom to use 
different representation of the BFKL cross section, equivalent within the NLO Q The good 
agreement between CMS data and BLM-optimized theoretical predictions gives us perhaps a 
hint towards the right direction. In the same paper, a list of issues was presented which, if 
considered in the experimental analysis, could help the matching between the way Mueller- 
Navelet are dehned in the theory and in the experiment. One of these issue was, for instance, 
the very measurement of the Mueller-Navelet total cross section, Cq, which, on the theory side, 
is strongly sensitive both to the representation of the BFKL amplitude and to the optimization 
procedure. 

In this paper we want to further discuss and expand another of the issues listed in the 
“Discussion” section of 123, related with the choice of the experimental cuts in the values of 
the forward jet transverse momenta. Since the Born contribution to the cross section Co is 
present only for back-to-back jets, its effect is maximized when symmetric cuts are used; on 
the contrary, in the case of asymmetric cuts, the Born term is suppressed and the effects of the 
additional undetected hard gluon radiation is enhanced, thus making more visible the BFKL 
resummation, in comparison to descriptions based on the hxed-order DGLAP approach, in all 
observables involving Cq. 

For this purpose, we compare predictions for several azimuthal correlations and their ratios 
obtained, on one side, by a hxed-order DGLAP calculation at the NLO and, on the other side, 
by BFKL resummation in the NLA. 

To avoid misunderstanding we note that in what follows our implementation of the NLO 
DGLAP calculation will be an approximate one. We just use here NLA BFKL expressions 
for the observables that are truncated to the O (af) order. In this way we take into account 
the leading power asymptotic of the exact NLO DGLAP prediction and neglect terms that are 
suppressed by the inverse powers of the energy of the parton-parton collisions. Such approach 
is legitimate in the region of large Y which we consider here. The exact implementation of 
NLO DGLAP for Mueller-Navelet jets is important, because it allows to understand better the 
region of applicability of our approach, but it requires more involved Monte Garlo calculations 
(some hrst results were reported recently in |30j). 

To single out the only effect of transverse momentum cuts, we consider just one represen¬ 
tation of the Mueller-Navelet cross section (the exponentiated one) and just one optimization 
procedure (the BLM one, in the two variants discussed in [ 31 ] )■ 

The paper is organized as follows: in the next section we give the kinematics and the basic 
formulae for the Mueller-Navelet jet process cross section; in section we present our results; 
hnally, in section we draw our conclusions. 

^We remark that the impact of this uncertainty on azimuthal correlations and their ratios is much larger 
than the one resulting from the adoption of different jet algorithms at the NLO. 
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2 Theoretical setup 


In this section we briefly recall the kinematics of the process and the main formulae, referring 
the reader to previous papers [22112Z] for the omitted details. 

The process under exam is the production of Mueller-Navelet jets [T] in proton-proton 
collisions 


p{pi)+p{p 2 ) ^ jet(/cjJ +iet{kj^)+X , (1) 

where the two jets are characterized by high transverse momenta, kj_^ ~ kj^ 3> ^qcd large 
separation in rapidity; pi and p 2 are taken as Sudakov vectors satisfying Pi = pi = ^ and 
2 {P 1 P 2 ) = s. 


In QCD collinear factorization the cross section of the process Q reads 
da 


dx j^dxj^d'^kj^d'^kj^ 


E 




f 

'0 Jo 




where the i,j indices specify the parton types (quarks q = u,d, s, c, b; antiquarks q = u,d, s, c, b; 
or gluon g), fi {x, qip) denotes the initial proton PDFs; Xi ^2 are the longitudinal fractions of the 
partons involved in the hard subprocess, while xjj j are the jet longitudinal fractions; pi? is the 
factorization scale; ddij {xiX 2 S,pf) is the partonic cross section for the production of jets and 
X 1 X 2 S = s is the squared center-of-mass energy of the parton-parton collision subprocess (see 

Fig.|5- 


The cross section of the process can be presented as 


da 

dypdyj^ d\kj^\ d\kj^\d(f)j^d(j)j^ 


1 

(27r)2 


00 

Co+ ^^2 cos(n0) Cn , 

n=l 


(3) 


where (j) = (l)j^ — (f)j^ — tt, while Co gives the total cross section and the other coefficients Cn 
determine the distribution of the azimuthal angle of the two jets. In the BFKL approach 
several NLA-equivalent expressions can be adopted for Cn- A large list of them and of their 
features can be found in [27] . For the purposes of the present analysis, we concentrate on one 
representation, the so-called exponentiated representation, and use it in combination with the 
BLM optimization procedure. We recall that BLM setting means choosing the scale pr such 
that it makes vanish completely the /do-dependence of a given observable. As discussed in [27] , 
we implement two variants of the BLM method, dubbed (a) and (6), derived in [21]. Moreover, 
we use a common optimal scale for the renormalization scale pr and for the factorization scale 
Pf- In [27] it was shown that this setup allows a nice agreement with CMS data for several 
azimuthal correlations and their ratios in the large Y regime. 


Introducing, for the sake of brevity, the definitions 


Y 


In 


^ Jl^ J2^ 

\kp\\kj^ 


Yo = \n ^ ^ 

\kjA\kj2 
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we have then the following expressions for the coefficients Cn-, in the two variants of BLM 
setting: 






\kjA\kj^ 


X {fin) ci{n,ii, |fcjJ,xjjc2(n, z/, \kj^\,xj^) 


( 4 ) 


X 


1 - as {fin) + as {fin) 


c?{n,iy, ^ 4"^(n, u, 

Ci(n, z/, C 2 (n, z/, \kj^\,xj^) 


.(i)z 


TT 


with fin hxed at the value 






( 5 ) 


and 


^BFKL(b) _ Xj^Xj^ 


^+oo 


^iy-yQ)[o!sifJ.R)xin,i')+a‘iilJ.R){x(n,i')-^xin,i^))] 


\kji I \ kj2 1 J 

X (fin) ci{n, u, \kn\, xj^)c2{n, u, \kj2\,xj2 


( 6 ) 


X 


/3q ^ 

1 + tts (hi?) ( {n, 1^) - 2 — j + as {fin) 


(n, z/, I fcjJ, X jJ ^ (n, z/, I A:j 2 1, X j 2 ) 


with /iji; hxed at the value 

(pBLM)2 ^ g^p 


ci(n, z/, |A;jJ,a;jJ C2(n, z/, |/cj2|,a;jJ 




In Eqs. 0 and 0, as{fin) = Cis{fin)Nc/'x^ with the number of colors, 


11 2 
/^o = y - ^nf 


is the hrst coefficient of the QCD /3-function, 


X (n, z/) = 2 z/> ( 1 ) - -0 ( ^ + ^ + w j + 


2 2 


( 7 ) 


( 8 ) 


(9) 


is the LO BFKL characteristic function. 


Ci{n,iy,\k\,x) = 2j^{Py‘' (^fg{x,fiF)+^fg{x,fiF)j (10) 

\ / 


and 


C 2 {n,i',\k\,x) = ci(n, z/, |fc|,x) 


( 11 ) 


are the LO jet vertices in the z/-representation. The remaining objects are related with the NLO 
corrections of the BFKL kernel {x{n-, z^)) and of the jet vertices in the small-cone approximation 
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|fcji 2 1, a:ji 2 )) in the z/-representation. Their expressions are given in Eqs. (23), (36) 
and (37) of Ref. Moreover, 


where/ = —2 Jq ~ 2.3439 andc[^ 2 (’^) \kj2\,xj^) are the same as 4^2(^) I^Ji, 2 1) ^Ji. 2 ) 

with the terms proportional to (3o removed. The scale sq entering ho is artihcial. It is introduced 
in the BFKL approach at the time to perform the Mellin transform from the s-space to the 
complex angular momentum plane and cancels in the full expression, up to terms beyond the 
NLA. In the following it will always be hxed at the “natural” value Tq = 0. 

In the hxed-order DGLAP approach at the NLO, the expressions for the coefficients Cn are 
nothing but the truncation of the BFKL expressions up to inclusions of NLO terms and read 

DGLAP(,) Xj^Xj^ J 2/ ^ / iT I ^ / iT I ^ 

\kjC\kj 2 \ /-CXD 


X 


1- as (hr) + as (hr) {Y - Yq) x {n, v) 


TT 


! c^l\n,v,\kjXxj^) c'2\n,v,\kj^\,xj^) 

+as {^^R) — -^-- + —-^-- 

\ ci(n, I/, C 2 [n,uCkj^\,xj^) 

XJ 1 XJ 2 f . 2 


Ai), 


0 


DGLAP 


(b) 


\kji\\kj2 


dv al {ij.r) Cl (n, I/, I /I, a; )C 2 (n, i/, | A;J 2 1, a; J 2 ) 


X 


1 + tts (hi?) ( {n, v) - 2— j + (hi?) {Y - Yo) X {n, v) 


(hi?) 


cS^^(n,z/, |/cji|,a;jJ ^ 4"^(n, z/, |A;j2|, a;j2) 


hi)/ 


( 12 ) 


(13) 


Ci(n, z/, |A;jJ,a:jJ C 2 (n, z/, |fcj 2 |,a:j 2 ) 

which we will use in the following with the two possible choices (a) and (6) of the optimal scales, 
given in Fqs. ([^ and ([^, respectively. It is worth mentioning that our DGLAP expressions, 
(12) and (13), do not actually depend on Yq. The corresponding terms in the r.h.s. of (12) 
and (13) are cancelled by similar terms in c^^ 2 ) 


We note that, in our way to implement the BLM procedure (see [31]), the hnal expressions 
are given in terms of as in the MS scheme, although in one intermediate step the MOM scheme 
was used. 
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Numerical analysis 


3.1 Results 

In this Section we present our results for the dependence on the rapidity separation between 
the detected jets, V = yj^ — yj 2 , of ratios TZnm = ^nl^m between the coefficients C„. Among 
them, the ratios of the form have a simple physical interpretation, being the azimuthal 
correlations (cos(m/))). 

In order to match the kinematic cuts used by the CMS collaboration, we will consider the 
integrated eoeffieients given by 

^yi,max ry2,iaax pOO pOO 

Cn= dyi / dy2 / dkj^ / dkjJ {yi -y2-Y) Cn {yj^.yj^, kj^,kj^) (14) 

^!yi,min ^?y 2 ,inin ^^Jj^,min ^J2,ra\n 

and their ratios Rnm = Cn/Cm- We will take jet rapidities in the range delimited by = 

1/2,min = —4.7 and 1 / 1 , max = l/2,max = 4.7 and consider F = 3, 6 and 9. As for the jet 
transverse momenta, differently from all previous analyses, we make two asymmetric choices: 
(1) fcji,min = 35 GeV, /cj2,min = 45 GeV and (2) = 35 GeV, kj^^^ain = 50 GeV. The jet 

cone size R entering the NLO-jet vertices is fixed at the value R = 0.5, the center-of-mass 
energy at ^/s = 7 TeV and, as anticipated, Yq = 0. We use the PDF set MSTW2008nlo [32] 
and the two-loop running coupling with as {Mz) = 0.11707. 

We summarize our results in Tables and and in Figs. [^and|^ We can clearly see that, 
at F = 9, BFKL and DGLAP, in both variants (a) and (6) of the BLM setting, give quite 
different predictions for the all considered ratios except C\/Cq] at F = 6 this happens in fewer 
cases, while at F = 3 BFKL and DGLAP cannot be distinguished with given uncertainties. 
This scenario is similar in the two choices of the transverse momentum cuts. 


3.2 Used tools 

All numerical calculations were implemented in Fortran, using the corresponding interfaces 
for the MSTW 2008 PDFs [32]. Numerical integrations and the computation of the polygamma 
functions were performed using specific CERN program libraries [33] . Furthermore, we used 
slightly modified versions of the Chyp [3l| and Psi [35| routines in order to perform the cal¬ 
culation of the Gauss hypergeometric function 2 F 1 and of the real part of the 'ip function, 
respectively. 

^In |27j it was mistakenly written 2/i,min = 0, although all numerical results presented there were obtained 
using the correct value for 
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3.3 Uncertainty estimation 

The are three main sources of uncertainty in our calculation; 


• The first source of uncertainty is the numerical 4-dimensional integration over the vari¬ 
ables l^jal, 2/ji and z/ and was directly estimated by Dadmul integration routine [HB] . 


• The second one is the one-dimensional integration over the longitudinal momentum frac¬ 
tion C entering the expression for the NLO impact factors l^ji 2 I) 2 ) given in 

Eqs. (36) and (37) of Ref. |22] and used in this work. This integration was performed 
by using the WGauss routine |33]. At hrst, we hxed the best value of the input accuracy 
parameter EPS by making comparisons between separate Fortran and Mathematic A 
calculations of the impact factor. Then, we verihed that, under variations by factors of 
10 or 1/10 of the EPS parameter, the and coefficients change by less than 

1 permille. 


• The third one is related with the upper cutoff in the integrations over |fcji|, {kj^l and u. 
We fixed fcjjmax = kj^miax = 60 GeV as in [21], where it was shown that the contribution 
to the integration from the omitted region is negligible. Concerning the //-integration, 
we fixed the upper cutoff z/max = 30 for the calculation of the coefficients, after 

verifying that a larger value does not change the result in appreciable way. 


The coefficients show a more pronounced sensitivity to z/max, due to the fact that 

the oscillations in the integrand in Eqs. (12) and (13) are not dumped by the exponential 
factor as in the BFKL expressions (|^ and For the same reason, the computational 

time of is much larger than for We found that the best compromise was 

to set z/max = 50. We checked in some sample cases, mostly at F = 6 and 9, that, putting 
z /max at 60, ratios Cm/Cn change always less than 1%, in spite of the fact that the single 
coefficients Cn change in a more pronounced way. 


Of the three main sources of uncertainty, the hrst one is, by far, the most signihcant, 
therefore the error bars of all data presented in this work are just those given by the Dadmul 
integration. We checked, however, using some trial functions which mimic the behavior of the 
true integrands involved in this work, that the error given by the Dadmul integration is a large 
overestimate of the true one. We are therefore conhdent that our error estimation is quite 
conservative. 


4 Conclusions 


In this paper we have considered the Mueller-Navelet jet production process at LHC at the 
center-of-mass energy of 7 TeV and have compared predictions for several azimuthal correlations 
and ratios between them, both in full NLA BFKL approach and in fixed-order NLO DGLAP. 
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Differently from current experimental analyses of the same process, we have used asymmetric 
cuts for the transverse momenta of the detected jets. In particular, taking one of the cuts at 
35 GeV (as done by the CMS collaboration [22]) and the other at 45 GeV or 50 GeV, we 
can clearly see that predictions from BFKL and DGLAP become separate for most azimuthal 
correlations and ratios between them, this effect being more and more visible as the rapidity 
gap between the jets, Y, increases. In other words, in this kinematics the additional undetected 
parton radiation between the jets which is present in the resummed BFKL series, in comparison 
to just one undetected parton allowed by the NLO DGLAP approach, makes its difference and 
leads to more azimuthal angle decorrelation between the jets, in full agreement with the original 
proposal of Mueller and Navelet. 

This result was not unexpected: the use of symmetric cuts for jet transverse momenta 
maximizes the contribution of the Born term, which is present for back-to-back jets only and 
is expected to be large, therefore making less visible the effect of the BFKL resummation. 
This phenomenon could be at the origin of the instabilities observed in the NLO hxed-order 
calculations of |3Sl 137] . 

Another important beneht from the use of asymmetric cuts, pointed out in [2S], is that 
the effect of violation of the energy-momentum conservation in the NLA is strongly suppressed 
with respect to what happens in the LLA. 

In view of all these considerations, we strongly suggest experimental collaborations to con¬ 
sider also asymmetric cuts in jet transverse momenta in all future analyses of Mueller-Navelet 
jet production process. 
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Table 1: Ratios Cn/Cm for /cji^min = 35 GeV and /i:j 2 ,min = 45 GeV. 



Y 

BFKL(,) 

DGLAP(,) 

BFKL(,) 

DGLAP(fe) 


3.0 

0.963(21) 

1.003(44) 

0.964(17) 

1.021(78) 

Gi/Go 

6.0 

0.7426(43) 

0.884(61) 

0.7433(30) 

0.914(91) 


9.0 

0.897(15) 

0.868(16) 

0.714(10) 

0.955(50) 


3.0 

0.80(2) 

0.948(43) 

0.812(15) 

0.949(75) 

G 2 /G 0 

6.0 

0.4588(32) 

0.726(56) 

0.4777(26) 

0.702(81) 


9.0 

0.4197(79) 

0.710(15) 

0.3627(50) 

0.850(48) 


3.0 

0.672(18) 

0.876(41) 

0.684(13) 

0.838(70) 

G 3 /G 0 

6.0 

0.3095(26) 

0.566(45) 

0.3282(21) 

0.435(68) 


9.0 

0.2275(72) 

0.558(13) 

0.2057(29) 

0.717(44) 


3.0 

0.831(18) 

0.945(43) 

0.842(16) 

0.929(72) 

G 2 /G 1 

6.0 

0.6178(43) 

0.821(66) 

0.6427(34) 

0.768(91) 


9.0 

0.4677(63) 

0.817(18) 

0.5079(56) 

0.890(51) 


3.0 

0.839(22) 

0.924(45) 

0.843(17) 

0.883(76) 

G 3 /G 2 

6.0 

0.6745(64) 

0.780(71) 

0.6869(52) 

0.62(11) 


9.0 

0.542(15) 

0.787(21) 

0.5670(59) 

0.844(56) 


Table 2: Ratios Cn/Cm for /cji,min = 35 GeV and 

min = 50 GeV 


Y 

BFKL(,) 

DGLAP(,) 

BFKL(,) 

DGLAP(fe) 


3.0 

0.961(23) 

1.006(46) 

0.964(15) 

1.034(89) 

Gi/Go 

6.0 

0.7360(49) 

0.869(58) 

0.7357(25) 

0.89(12) 


9.0 

1.0109(61) 

0.857(16) 

0.7406(46) 

0.958(56) 


3.0 

0.788(21) 

0.946(44) 

0.801(14) 

0.950(85) 

Gs/Go 

6.0 

0.4436(37) 

0.698(53) 

0.4626(19) 

0.611(98) 


9.0 

0.4568(50) 

0.695(15) 

0.3629(23) 

0.862(54) 


3.0 

0.653(19) 

0.868(43) 

0.669(12) 

0.814(79) 

G3/G0 

6.0 

0.2925(31) 

0.530(42) 

0.3115(15) 

0.320(57) 


9.0 

0.2351(35) 

0.551(17) 

0.1969(17) 

0.748(50) 


3.0 

0.820(21) 

0.940(44) 

0.832(15) 

0.918(81) 

G2/G1 

6.0 

0.6027(51) 

0.803(64) 

0.6288(26) 

0.69(12) 


9.0 

0.4518(35) 

0.811(18) 

0.4900(24) 

0.899(57) 


3.0 

0.829(26) 

0.917(46) 

0.835(17) 

0.857(85) 

G3/G2 

6.0 

0.6595(82) 

0.759(70) 

0.6733(36) 

0.52(11) 


9.0 

0.5146(85) 

0.793(23) 

0.5426(38) 

0.869(62) 
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Figure 2: y-dependence of several ratios Cm/Cn for fcj^,min = 35 GeV and fcj 2 ,min = 45 GeV, 
for BFKL and DGLAP in the two variants of the BLM method (data points have been slightly 
shifted along the horizontal axis for the sake of readability). 
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Figure 3: y-dependence of several ratios Cm/Cn for fcj^,min = 35 GeV and = 50 GeV, 

for BFKL and DGLAP in the two variants of the BLM method (data points have been slightly 
shifted along the horizontal axis for the sake of readability). 
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